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Abstract 

We derive a new parametric class of exact cosmological solutions to Brans-Dicke theory of 
gravity with a self-interacting scalar field and a barotropic perfect fluid of ordinary matter, 
by assuming a linear relationship between the Hubble expansion parameter and the time 
derivative of the scalar field. As a consequence only a class of exponential potentials and 
their combinations can be treated. The relevance of the solutions found for the description of 
the cosmic evolution are discussed in detail. We focus our discussion mainly on the possibility 
to have superquintessence behavior. 



1 Introduction 



Dark energy or missing energy is one of the contemporary issues the physics community 
is more interested in due, mainly, to a relatively recent (revolutionary) discovery that our 
present universe is in a stage of accelerated expansion[l]. This missing component of the ma- 
terial content of the universe is the responsible for the current stage of accelerated expansion 
and accounts for 2/3 of the total energy content of the universe, determining its destiny[2]. 
This is a new form of energy with repulsive gravity and possible implications for quantum 
theory and supersymmetry breaking. Among others, the new cosmology is characterized 
by the following features[2]: 1) Flat, critical density accelerating universe, 2) early period 
of inflation, 3) density inhomogcncitics produced from quantum fluctuations during infla- 
tion, 4) composition of 2/3 dark energy, 1/3 dark matter and 1/200 bright stars, 5) matter 
content: (29 ± 4)% cold dark matter, (4 ± 1)% baryons and ~ 0.3% neutrinos. Besides, 
there is evidence (recent observation of SN1997ff at redshift z = 1.7) that the present stage 
of accelerated expansion was preceded by an early period of decelerated expansion[2]. Any 
theoretical framework intended to describe our present universe should be compatible with 
these observational evidences. 

A self-interacting, slowly varying scalar field, most often called quintessence, has been 
meant to account for the dark energy component. In the simplest models of quintessence, 
minimally coupled scalar fields are considered. However, extended quintessence models have 
also been considered[3, 4, 5]. In these models the scalar field that accounts for the dark 
energy component is non-minimally coupled to gravity and, besides, is a part of the gravity 
sector. Brans-Dicke theory is a prototype for this kind. Motivated by the possibility to 
have superaccelerated expansion and, correspondingly, superquintessence solutions, this type 
of theories is now being considered (see, for instance, Ref.[5] and references therein). In 
particular, in Ref . [5] . the author considered different non-minimally coupled theories (among 
them the Brans-Dicke theory) and showed that, although superquintessence may arise in 
these models, its effect is vanishingly small, i. e., is far bcyong the present observational 
possibilities. In that paper the author numerically integrated the field equations by using a 
computed code. Instead, it would be of interest to have a parametric class of exact solutions 
to treat this subject. 

The aim of the present paper is, precisely, to derive a new parametric class of exact 
cosmological solutions to Brans-Dicke gravity theory with a self-interacting scalar field, by 
assuming a linear relationship between the Hubble expansion parameter and the time deriva- 
tive of the scalar field. This type of relationship is a particular case of a more general one 
that has been used in Ref. [6] to find general solutions in N + 1-dimensional theories of 
gravity with a self-interacting (minimally coupled) scalar field and also to derive general 
solutions when we have a self-interacting scalar field with exponential potential plus a free 
scalar field. ^ Also, it has been already used by us[8] to derive exact cosmological solutions 
to Einstein's gravity with two fluids: a barotropic perfect fluid of ordinary matter, together 

method where a linear relationship between the field variables and/or their derivatives is assumed has 
also been used in [7] to derive 4d Poincare invariant solutions in thick brane contexts. 
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with a self-interacting scalar field fluid accounting for the dark energy in the universe. The 
assumed relationship between the Hubble parameter and the time derivative of the scalar 
field is suggested by an implicit symmetry of the field equations. 

We point out that it is not neccessary to make any a priori assumptions about the func- 
tional form of the self-interaction potential or about the scale factor behavior. These are 
obtained as outputs of the assumed linear relationship between the Hubble expansion pa- 
rameter and the time derivative of the scalar field, once one integrates the field equations 
explicitely. As a consequence only a class of exponential potentials and their combinations 
can be treated. However, this is not a serious drawback of the method. In effect, combi- 
nations of exponentials are interesting alternatives since these arise in more fundamental 
(particle) contexts: supergravity and superstring[9], where these types of potentials appear 
after dimensional reduction. The relevance of the solutions found for the description of 
the cosmic evolution will be discussed in some detail making emphasis in the possibility 
for superquintessence solutions to arise. This issue has been alredy studied numerically in 
Ref.[5]. 

We will be concerned with flat Priedmann-Robertson- Walker (FRW) cosmologies with 
the line element given by: 

ds"^ = -dt^ + a{tf6ikdx'dx'', (1) 

where the indexes i,k = 1, 2, 3 and a{t) is the scale factor. We use the system of units in 
which SttG = c = 1. 

2 The Method 

The Brans-Dicke field equations, in the presence of a self-interacting scalar field, are: 

3H^ + 3Hip - |</.2 = e-'^ipm + V), (2) 

2H + 3H^ + ^ + ip^ + 2Hip + ^ _ + V, (3) 

+ <p' + 3H<p = ^^{(4 - 3l)Pm + 4y - 2V'}, (4) 

where we have introduced a new field variable (p instead of the original (p — e'^. In these 
equations cu is the Brans-Dicke coupling parameter, 7 is the barotropic index of the fluid 
of ordinary matter, V is the self-interaction potential and H = - is the Hubble expansion 
parameter. The dot accounts for derivative in respect to the cosmic time t meanwhile 
the comma denotes derivative in respect to the scalar field ip. The energy density of the 
ordinary matter (cold dark matter plus baryons and/or radiation) is related with the scale 
factor through — po^^a~^'^, where po,7 is an integration constant coming from integrating 
the conservation equation. 
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Let us introduce a new time variable dt — e'^/^dr. Then the field equations (2-4) can be 
rewritten in the following way: 



3H^ = pm + Pv', (5) 
2H + = -Pm - P^, (6) 



where 



Pm = {l- l)Pm. (8) 

Now the dot means derivative in respect to the new time variable r and H — 
Combining of the above equations yields 



fi + % + 3i?'^=<i^2^p„ + 5^y + ^y'. (9) 

2 ^ 2u; + 3 2a; + 3 2a; + 3 ^ ^ 

An implicit symmetry of the LHS of equations (7) and (9) is made explicit by the change 
H kif, where k is a constant parameter. Therefore, if we assume the linear relationship 

H = kip, a = e'"^, (10) 
then, the following differential equation for the self-interaction potential V is obtained: 

2a; + l-4A; (4 - 37)A: - (2 - 7)^ - 1 ^ .... 
^ 2A; + 1 2A; + 1 ^'"^ ^ ^ 

Straightforward integration of Eq. (11) yields 

V{ip{a))^Aa'' + Ba-^\ (12) 

where 

4A; - 2a; - 1 

n - 



{2k + l)k ' 



37A; — a 6^k — a 

, . (4-37)fc-(2-7)a;-l 
a^-kn, /3 = . (13) 
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An interesting feature of this potential is that it depends on the type of ordinary fluid 
which fills the universe. Otherwise, it depends on the barotropic index 7 of the matter fluid. 
This fact implies some kind of interaction between the ordinary matter and the quintessence 
field much like the interacting quintessence studied in Ref.[10]. 

By substituting Eq. (12) back into the Priedmann equation (5) one obtains 

= Aa" + (14) 

where 

6F + 6A;-a;' + 6A; - a; ' ^ ^ 

Finally we make another change of time variable: dr — oP'/'^dT to get Eq. (14) integrated 
in quadratures: 

3-T+n 1 , 

= r + ro, (16) 



where Tq is another integration constant. 

3 The Class of Solutions 

Straightforward integration of Eq. (16) leads to: 

a(r) = ao{sinh[/i(r + ro)]}^^, (17) 

with 

ao = (^/i)^, = (18) 
Then, from the linear relationship (10), one can derive the form of the scalar field: 

1 2 

(^(r) = -lnao+ ^ lnsinh[//(r + tq)]. (19) 

For purpose of observational testing we give the list of magnitudes of observational in- 
terest. The Hubble expansion parameter; 

H{r) = ^ a{rY'^ coth[//(r + ro)], (20) 
07 + n 

the matter density parameter, 

^m{r) = ^ = (^)'^ cosh-[/.(r + ro)], (21) 



the scalar field density parameter. 
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Q^(r) = l-Q^(r), (22) 

the deceleration parameter, 

?W = -l + ^^ + ^Mr), (23) 

and the standard (perhaps not precisely the physically meaningful [5]) equation of state of 
the scalar field, 

u^^^^-l- n/3. (24) 
This class of solutions depends (in principle) on six parameters: 7, cu, k, po,7, C'o and tq. 



4 Reducing the Space of Parameters 

A very crude analysis shows that should never be negative since, otherwise, this would 
imply coarse violations of the known energy conditions. Therefore should never be 
greater than unity Since flm (Eq- (21)) is a maximum at r = — ro where 



iC(-ro) = ( ) (25) 



then by requiring that 



M-ro) = l, ^ ^ = J— • (26) 
If we apply the normalization in which po,7 — 1) then = ^/3 and Eq. (26) implies: 

_ (67 + (4 + 37)/c - 2cj - l){Qk^ + 6k-uj) 
^° " 6A;2(2A: + 1) " 

By considering Eq. (26) and evaluating Eq. (21) at present (r = 0) one obtains: 

To = —arccoshl , =], (28) 



where, from observations, one may place Qm{0) = 0.3. 

After these considerations the space of parameters is reduced to the following free pa- 
rameters: 7, uj and k. Besides the parameter 7 is fixed by the type of fiuid dominating the 
given stage of the cosmic evolution. At present (r = 0) 7 = 1, i. e., the universe is filled 
with a cold dust of dark matter plus baryons mainly. 

We want to point out that all of the physical parameters of observational interest (equa- 
tions (21-24)) can be given as functions of the redshift z simply by replacing in the former 
equations 



5 



This form of writing will be useful for purpose of observational testing. 



5 Observational Testing 

The main observational facts we consider are the following [2]: 

1. - At present (r = 0) the expansion is accelerated (g(0) < 0). 

2. - The accelerated expansion is a relatively recent phenomenon. Observations point to 
a decelerated phase of the cosmic evolution at redshift z = 1.7. There is agreement in that 
transition from decelerated into accelerated expansion occurred at a 2; ~ 0.5[11]. 

3. - The equation of state for the scalar field at present a;,^ ~ — 1 (it behaves hke a 
cosmological constant). 

4. - Although, at present, both the scalar (quintessence) field and the ordinary matter have 
similar contributions in the energy content of the universe ($1^(0) = 0.3 ^ f2ip(0) = 0.7), 
in the past, the ordinary matter dominated the cosmic evolution,^ meanwhile, in the future, 
the quintessence field will dominate (it already dominates) and will, consequently, determine 
the destiny of the cosmic evolution. 

In fig.l the evolution of both the matter and scalar field density parameters is shown 
for ijj = 1000 and an arbitrary k = 530 (7 = 1), as function of the redshift z. We see that 
the evolution of both magnitudes meets the observational requirements o matter dominance 
(scalar field subdominance) at high redshift (in the past) and scalar field dominance at low 
redshift (at present and in the future). 

Prom fig.2, where the evolution of the deceleration parameter g & 2; is shown for the 
same u = 1000 and k = 530, we see that, at present {z = 0), g(0) = —0.55 < (accelerated 
expansion) meanwhile, in the past {z > 0) we had a stage of decelerated expansion. The 
transition from decelerated into accelerated expansion occurred at z = 0.67, which is in 
acceptable agreement with the claimed value z 0.5. Besides, the observational fact of 
decelerated expansion at redshift z — 1.7 is fulfiled. 

The equation of state is always a constant u^p = —1.00007 for the former values of the 
free parameters uj and k, in good agreement with the accepted value. In general the solution 
found meets all of the above observational requirements (points 1-4) for a really wide range 
in the parameter space {u!,k). 



6 Constrains from Superquintessence 

As pointed out in Ref. [3], superquintessence behavior {u^ < —1) is a feature of models of 
quintessence based on non-minimally coupled scalar-tensor theories including Brans-Dicke 

^The quintessence field (fi sliould be subdominant in the past as required by nucleosynthesis constraints[12] 
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theory with a self- interacting scalar field. In Ref.[5], however, it was shown that expected 
superacceleration effect is small enough as to escape the reach of present and future obser- 
vations. We now proceed to analyze the question of superquintessence within the frame of 
the solution found. 

We shall interested in n > (Eq. (13)) in order to have superquintessence behavior (see 
Eq. (24)). This, in turn, imphes k > (2a; -|- l)/4. Since, from experiment u; > 500[13], then 
k > 250.25 is a necessary condition to have superquintessence. Consider n as function of k 
{cu constrained by experiment): 

4:k-2uj-l 



which is a maximum n* at A; = (2a; -|- 1 -|- {2cu + l){2u) + 3))/4, where 



n* = 2V(2c. + l)(2.; + 3) 

{^(1 + ySf )P + (^)(i + ySf ) ■ 

For instance, for lo — 500, the maximum value for n is n* ^ 0.002 (at k ^ 500.75) so, 

according to Eq. (24), the maximum expected superacceleration effect is of very small magni- 
tude (0.2%). For other values of /c > {2uj+ l)/4, this effect would be even smaller. However, 
as discussed in details in Ref. [5] , the physically meaningful definition of the equation of state 
for the scalar field reads; 

l ln[p^(r)/p^(0)] 
^^-"^"3 ]^) ' 

instead of Eq. (24). In (32) the following redefinition for the scalar field energy density has 
been used: 

p^^p^- SH^{1 - oo '/^(sinh[/i(r + ro)])~(OT}. (33) 

In fig. 3 we show the behavior of the equation of state both, the usual (Eq. (24)) and the 
physically meaningful (Eq. (32)) [5], as function of the parameter A; for a; = 1000. Maximum 

expected superquintessence, in this case, is obtained for k ~ 800 (the standard equation of 
state) or A; ~ 1300 (the physically meaningful equation of state). In both cases the effect is 
very small (of about 0.( 



7 Conclusions 

We have found a new parametric class of exact cosmological solutions to Brans-Dicke theory 
of gravity with a self-interacting scalar field. A class of self-interaction potentials that is a 
combination of exponentials and that is dependent on the kind of matter fluid that fllls the 
space-time arise. This potential suggests some kind of interaction between the quintessence 
field and the ordinary matter similar to the kind in Ref. [10]. 
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The flat universe solution found meets the observational constrains on the physical pa- 
rameters, among them Q^n (^ip)) ^(^j Q ^-^d the experimental requirement that u > 500 [13]. 
It reflects the observational fact of a past decelerated, ordinary matter (including radiation) 
dominated universe and a future scalar held (dark energy) dominated universe. 

The analytic study of the class of solutions derived in the present paper conflrm the results 
of Ref. [5]. Superquintessence behavior, although present for a given region in parameter 
space (oUjk), is behind the reach of present and future observational possibilities. 
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Figure 1: The matter density parameter Qrn (solid line) and the scalar field density parameter 
flfj) (dotted line) as functions of the redshift z for j = 1, uj = 1000 and an arbitrary k = 530. 
We have used the normalization in which po,7 = 1- We have considered the observational 
fact that, at present, 0^(0) = 0.3. It is seen an early stage when the contribution from the 
quintessence field was subdominant (high redshift). At present {z = 0) both contributions 
from dust dark matter (plus baryons mainly) and from the scalar field are of the same order 
{flm{0) = 1/3 while ^2,^(0) = 2/3). In the future the quintessence field will be dominant. 



Figure 2: The deceleration parameter q is shown as function of the redshift z. At present 
qo = —0.55. It is seen that there are both an early decelerated phase where ^ > and a 
late accelerated expansion period (^ < 0). The transition from decelerated into accelerated 
expansion occurs at 2; = 0.67. 



Figure 3: The standard equation of state for the scalar field c<j<^ and the physically meaningful 
scalar field state equation uj^ are shown as functions of the free parameter k. We have chosen 
7 = 1 and uu = 1000 for definiteness. It is seen that both curves attain a minimum at some 
k, i. e., a maximum superquintessence effect takes place at these A;- values. This effect is very 
small and falls outside of the reach of present (and future) observational posibihties[5]. 
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